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Abstract 

We consider two classes of probability distributions for configura- 
tions of the ideal gas. They depend only on kinetic energy and they 
remain of the same form when degrees of freedom are integrated out. 
The relation with equilibrium distributions of Tsallis' thermostatistics 
is discussed. 



1 Introduction 

This paper studies two classes of probability distributions for the ideal gas 
with n degrees of freedom. Their main property is that elimination of some 
of the degrees of freedom produces again a distribution of the same class. 

Particle velocities in an ideal gas with infinitely many degrees of freedom, 
in thermodynamic equilibrium at inverse temperature j3, follow the Maxwell 
distribution. This result is known for more than a century. But quite often 
one is interested in finite systems or in systems out of equilibrium. For 
such systems deviations from a Gaussian distribution are expected. Such 
deviations have been reported recently pQ in laser-cooled dilute samples of 
atomic gases. The velocity distribution of migrating Hydra cells in two- 
dimensional aggregates is non-Gaussian, the mean-square displacement is 
superdiffusive j2]. Numerical simulation of metastable states in a long-ranged 
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mean-field model j3 IH El reveal non-Gaussian velocity distributions. Many 
more examples should be known but a systematic investigation seems to be 
missing. 

If the velocity distribution is not Maxwellian then one cannot expect it 
to be of the product form. Indeed, the standard argument 



works only for the exponential function. This raises the question whether 
one can formulate alternatives which are still easy to manage. Distributions 
(JTEJUBJ), studied in the present work, satisfy the requirement of simplicity. 
They depend only on the total energy of the gas. They have the additional 
property that, if some of the degrees of freedom are integrated out, then 
the remaining degrees of freedom obey a probability distribution of the same 
form. Two classes of distributions are introduced. They correspond with the 
cases q > 1 and q < 1 in the terminology of non-extensive thermostatistics. 
The first class contains distributions for which the probability of large ener- 
gies decays algebraically instead of exponentially. The other class consists of 
distributions with a cutoff for large energies. 

In the next section the two classes of probability distributions are de- 
rived starting from simple models for the velocity variables. In Section 3 
is shown that the distributions depend only on total energy and a duality 
relation between the two classes is pointed out. In Section 4 the relation with 
equilibrium distributions of non-extensive thermostatistics is explained. The 
latter do not generically satisfy the requirement of stability under reduction 
of degrees of freedom. In the final section follows a discussion of our results. 

2 Two families of probability distributions 

A characteristic feature of the ideal gas with n degrees of freedom is that 
it is described by independent variables k\, k 2 , ■ ■ ■ , k n , components of the 
velocities of the particles. It is tradition to assume that each of the variables 
kj obeys the Maxwell distribution 



with (3 the inverse temperature and with rrij the mass of the j-ih degree 
of freedom. A slight generalization is obtained by replacing the Gaussian 
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distribution (X/tt) 1 ^ 2 exp(— Xu 2 ) by an exponential power distribution 

f^(«) = 2^^«p(-A|tt|*) ) (3) 

with A the scale parameter and z the shape parameter. These distributions 
were introduced by Subbotin [S| in 1923. The distribution (J2J) is then replaced 

by V^%{kj) with 

Xj = -£m,, (4) 

This generalization is not essential for the paper. One can take z = 2 in all 
what follows. 

Several arguments can be invoked to explain why sometimes the distri- 
bution of experimentally observed velocities Vj differs from that of the in- 
dependent variables kj. The simplest argument is that the distribution is 
modulated by an independent variable s. The relation between v and k is 
assumed to be of the form 

v i = s x /z ( 5 ) 

with g{0) a positive proportionality constant, depending on inverse temper- 
ature f3. We will assume that s obeys a Gamma-distribution with shape 
parameter 6+1 



The probability distribution of the velocities v is then found to be 

p+oo n 

Pn(v)= dsw(s) \\pl:%/g(p){vj). (7) 

J ° 3=1 

An interesting alternative for (J3J) is 

V , = (8) 



where s n has a Gamma-distribution with shape parameter b n . These proba- 
bility distributions form a family stable under reduction of degrees of freedom 
because 

s n -i = s n + X n \k n \ z (9) 
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has again a Gamma-distribution, this time with shape parameter b n _i given 
by 

b n -i = b n + ~ (10) 

For further use let b n = b — l — n/z. The resulting probability distribution for 
the velocities Vj, given by (JEJ), is denoted q n {v) and, after some calculation, 
is found to be 



Qn{V) 



r+oo n 

/ ds n w n {s n )J n {s n ,k)T\p^ p x (k j ) (11 
Jo j=1 



with k related to v by (jBJ), with 



w - w = fh st " le " (12) 



and with 



J.(*,*)=det(*W , . (13) 

3 Energy distribution and duality 

The energy of the ideal gas of n particles is given by 

^ n v 

E n{v) = - S ^m^\v ha \ z . (14) 

j=l a=l 

The probability distributions ((Zj) and fTTj) depend only on the velocities v 
via the energy E = H n {v). Indeed, a short calculation gives 

p " ( "»" p " (£) = ^m W)+ l m M ^' (15) 

q n (v)=q n (E) = ^ [<,(/?) - (16) 



with 
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UP) 



l-n/z) f2T(l/z)\ n 



T(b-1) 



b-1 



(17) 



The notation [u] + = max{0,ii} is used. Note that n < z(b — 1) is required 
for q n (v) to be well-defined. 

Let us now clarify the relation between p n {E) and q n (E). First note 
that expression (fT5|) is meaningful for negative n, as long as z(b + 1) + 
n > is satisfied. Similarly, (fTo^) is meaningfully for all negative n. The 
transformation 



E' 



g(P)E 



g(/3)+(3E 

maps p- n (E) onto q n (E'). Indeed, one has 



P-n{E) 



dE 

^P-n(E) 

g(P) 2 



(g((3) - (3E>y Z_ n ((3) 
1 (g(g) - 

Z_ n (/?) g([3f(b+l-n/ Z ) 

qn(E') 



g(P) + 



l—n/z 



Pg{P)E' 
g(P)-PE'\ 



-b-l+n/z 



(19) 



provided that the normalization constants Z n (J5) and ( n {P) for negative val- 
ues of n are defined in an appropriate way. 



4 Non-extensive thermostatistics 



In the context of Tsallis' non-extensive thermostatistics various proposals 
have been made for the equilibrium probability distribution at inverse tem- 
perature p. One of these reads [Ij 



p(E) 



1 - (1 - q)P 



Z(P) [ 

with normalization factor Z(P) and with 

1 



E - U q (P) 

z(py-« 
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u q (P) 



z(py- 



dEp(E) q E. 



(20) 



(21) 



In the limit that the parameter q equals 1 then (f2T)j) reduces to the well-known 
Boltzmann-Gibbs distribution. Identification of (f2T)|) with (|THj) gives 



g(P) 



z{P) 



q-l 



pu q (p) 
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Identification with (frfljl gives 



l-q 



1 - 
1 



+ pu q (p) 




l-q 



(23) 



These relations show that the parameter q must depend on the number of 
particles n. In general this identification is only possible for a single value of 
n because the function g(0) should not depend on n. 

Recently, one of the authors [BUS] has proposed equilibrium distributions 
of the form 



where <p(v) is an arbitrary positive function. Identification with (|151 ITfij) is 
only possible with a function <p which depends on (3 and on n. 

One concludes that in general the probability distributions used in non- 
extensive thermostatistics are not compatible with the reduction requirement 
which is central to the present work. 

5 Discussion 

Two classes of probability distributions for the velocities of an ideal gas have 
been introduced. The distributions p n (v) given by (fTHjl have the property 
that integrating out one of the velocities produces a probability distribution 
of the same class 



The same property holds for the distributions q n (v) given by (fTB|). In ad- 
dition, both kinds of probability distributions depend only on total kinetic 
energy. 




(24) 




(25) 




(26) 
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For large values of energy E the probability distribution p n (v) behaves as 



with constants b and z (usually is z = 2). Because of the n- dependence 
of the exponent the distribution of large systems approaches the Maxwell 
distribution. On the other hand, the probability distribution q n (v) vanishes 
for E larger than the cutoff value g((3)/(3. In addition, q n (v) is only defined 
for a number of particles less than some maximum value z(b — 1). Hence it 
is not possible to consider the thermodynamic limit n — > +00. 

Our derivation of the probability distribution p n (v) in Section 2 starts 
with a model of velocities modulated by a variable which obeys a Gamma- 
distribution — see expression (J3J). This same model is the basis of super- 
statistics which explained distributions of Tsallis' nonextensive thermo- 
statistics in case q > 1. For the case q < 1 no such explanation was available. 
Here we show how a modification of (0) leads to the distribution q n (v), which 
in the Tsallis context corresponds with q < 1. 

In Section 4 is shown that in the generic case the distributions of Tsallis' 
nonextensive thermostatistics do not satisfy the requirement (}26|) of stabil- 
ity under reduction. Note that in the context of finite-size renormalization 
a condition weaker than (|2T)|) is considered. This point of view has been 
elaborated in [TT] . 
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